Role of range of interactions in a model of diffusing charged particles 



Krzysztof Pawlikowski' EI and Andrzej P§kalski^'[l| 

'institute of Theoretical Physics, University of Wroclaw, 
pi. M. Borna 9, 50-204 Wroclaw, Poland 

Role of range of interactions in a model of charged particles diffusing on a two-dimensional lattice is studied. 
We investigate, via Monte Carlo simulations, three models. In the first one interactions are restricted to nearest 
neighbors, in the second one we add (weaker) interactions to second nearest neighbors. In the third model 
interactions are extended as far as the half the size of the system. The strength of these interactions falls off 
with the distance. We show that for static properties, such as asymptotic spatial patterns or the order parameter, 
the results of the first and third models are very close. Dynamic characteristics, like mean square displacement, 
average cluster size and average energy as functions of time show larger differences among the three models. 
Origins of these differences are pointed out and discussed, 
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I. INTRODUCTION 

Motion of particles and pattern formation on planar surfaces is a much studied subject fll 01 • Understanding basic mechanisms 
of such processes was always important, as it describes motion of adatoms or ions on surfaces. It became recently even more so 
with the growing interest in nanostructures IH]. Diffusion of particles on a substrate is a cooperative phenomenon resulting from 
interactions among many particles, particles with a substrate and from thermal excitations. Theoretical papers devoted to these 
problems are many. Some of them are based on analytic approach iB, 3, Bl , others used Monte Carlo simulations [7l[8il9l [Toll . 
The range of processes studied includes pattern formation IIQl lTolflltl or dynamics iQ-SSLIu some of these papers the studied 
systems was composed of one type of particles, in others there were two types, A and ijflSSl- Sometimes also vacant sites 
were allowed. In general, interactions of diffusing particles with the substrate atoms were neglected and the role of the substrate 
was reduced to creating an unchanging periodic potential acting on the diffusing particles. Forces between theparticles were 
modelled in many ways. Simple nearest-neighbor (NN) interactions fs*], next nearest-neighbors (NNN) ones fjj^ and even 
farther ||^[l3l- The interactions could be of one type (repulsive or attractive) or they could be of both types, depending on the 
distance between the interacting particles |8]. 

In two recent papers 1 12, 13] we have introduced a simple model of two types of particles moving on a two-dimensional square 
lattice. The particles are identical, except for having opposite charge. There is the same number of positively and negatively 
charged particles, hence the net charge is always zero and the dynamics is of the Kawasaki type. A particle is free to move to 
an empty NN site according to the Metropolis algorithm. Restricting interactions to NN and weaker NNN ones, we have shown 
how the resulting patterns and dynamic parameters, like diffusion coefficient, mean square displacement, average cluster size 
depend on the range of interactions (to NN or to NNN), coverage and the temperature. 

In this paper we investigate, using Monte Carlo simulations, a problem which is quite important for modeling diffusion on 
surfaces and which, to the best of out knowledge, has not yet been investigated, namely how important is the restriction of 
the interactions to NN and NNN. It is well known that if the local charge is screened, the range of the potential between 
interacting particles could be reduced. We shall compare here spatial patterns, order parameters, mean square displacement, 
tracer diffusion coefficient, average cluster size and average energy for three versions of the basic model introduced in lfl2ll . 
In the first version, which we shall henceforth call Hi, the interactions are restriced to NN, in the second {H2) one (weaker) 
interactions to NNN are added. In the third version, called H^, each particle interacts with all particles which are not farther 
from it than half the size of the system. Therefore, with periodic boundary conditions used by us, each particle interacts with 
all others, although the strength of such interactions falls off with the distance between the particles, contrary to mean-field 
models, where the interactions to all particles have the same strength. Quite often restricting the range of interactions of charged 
particles, acting with Coulomb forces, is criticized. It seems therefore useful to check how different would be the results coming 
from interactions reduced to NN, then to NNN, as compared with the long range ones. 
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II. MODEL 



A system of positively ( A) and negatively (B) charged particles located at sites of a square lattice of dimensions L x L is 
considered. The lattice has periodic boundary conditions and odd number of nodes. The total number of particles of both types 
divided by the number of sites defines the coverage, which in the following will be equal 0.2. A lattice site could be either empty, 
or occupied by one Aoi B particle. Hamiltonians for the three models described above have the form 

N Si 

Hi = Ji^^niUj , (1) 

where N is the total number of particles, Ji is a positive constant and rii = 0, ±1, depending whether the site i is empty, 
occupied by an A or B particle. 6i means summation over first neighbors. When the interactions are extended to NNN the 
Hamiltonian takes the form 

N / Si S2 \ 

^2 = ^ "-j Ml ^ nj + ^2 ^ , (2) 

i=l \ J = l J = l / 

where J2 = Ji/V^ and S2 means summation over NNNs. Finally, when the interactions are up to any particle within the 
{L — l)/2 range, we have 

JV 

7^3=^ niUj. (3) 

Now the sum extends to neighbors not farther than (i — l)/2 and the interaction constant is divided by the distance between the 
sites i and j. The three models belong to the class of stochastic lattice gas models with additional degree of freedom (charge). 
In one Monte Carlo step (MCS), which is our time unit, each particle could move to an empty NN site according to the standard 
Metropolis algorithm jlSll . Constant temperature in the system is maintained through a contact with a heat reservoir 

We shall investigate the following quantities: the order parameter, asymptotic, i.e. at the end of simulations, spatial pat- 
terns, mean square displacement (msd), average energy and average cluster size as functions of time. From the mean square 
displacement {r^{t)) the tracer diffusion coefficient D* is calculated as llj] 



where fi (t) and (0) are the positions of the particle i at time t and time 0, respectively. Another quantity recorded is the order 
parameter (p introduced in Ref. lfT2ll and generalized to interactions up to m neighbors as 
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i=i i=i 



(5) 



This order parameter is zero in high temperatures (no stable pattern) and is non-zero in low temperatures. Time dependence of 
the average cluster size yields useful information about the emerging spatial structure of the system and its stability. Following 
Landau and Binder [ 15] a cluster is an agglomeration of minimum 10 particles which have at least one particle as a NN or NNN. 
Average energy is the sum all interactions in the system divided by the twice the number of particles. 

As we have explained before ifTsll . we have studied a system of rather modest size, namely 51 x 51. In many cases to see 
interesting results we had to run simulations to quite long times, up to 3 x 10^ MCS. In large systems processes are going more 
slowly, hence to observe such effects one would have to run simulations for at least two orders of magnitude longer times, which 
would make it out of our computing powers. We have checked before, Ref. ifisll . that for low coveraged reported in this paper, 
there are no differences between systems with L = 51 and L = 101. Most of the results showing time dependence of average 
energy etc. come from single simulation runs. The reason is that characteristics of these curves are, in general, shifted in time 
in different runs. Averaging will result in smearing out of these effects. We have verified that the presented curves are generic 
and not just particular effects. All results presented below are for one, rather low, coverage c = 0.2. Discussion of the effects of 
higher coverages for the models Hi and H2 has been given in Ref. lfT2l[T3ll . Data for quantities like average energy, cluster size 
or msd are reported for three sets of temperatures - low (below the critical temperatures Tc), medium (a little above Tc) and high 
(much above Tc). The Arrhenius and order parameter plots are given for all interesting temperatures. 
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FIG. 1 : Order parameter versus scaled temperature for three models - interactions to NN, up to NNN and the long range ones. Scatter of the 
data is marked by error bars. 



FIG. 2: Spatial structures obtained at low temperatures after 4 x 10 MCS for the three models: Hi (left upper panel), H2 - (right upper 
panel), H3 (left bottom panel) and H3 after 3 x 10^ MCS (right bottom panel). 



III. RESULTS 



We begin with the order parameter ip. Quite often, especially in investigating critical phenomena, relevant parameter is not 
just the temperature, but rather its deviation from the critical value. Since in this paper we compare results of models differing 
by interactions an appropriate parameter will be the reduced (scaled) temperature, r — T/Tc, i.e. temperature divided by its 
critical value. From the dependence of on T we have estimated Tc as the inflection point of the (p{T) curve. The plot of (p{t) 
is shown in Figure [T] 

As could be seen, the scaling is quite good, showing the same, or at least very similar, mechanisms, regardless of the range of 
interactions. In the following we shall refer to temperatures r w 0.80 as low temperatures, to r w 1.02 as medium ones and when 
T « 1.60 as high temperatures. In absolute values Tc/Ji ~ 0.44, 0.31 and 0.081 for the Hi, H2 and models, respectively. 
From Fig[T]we see that a better agreement is reached between the Hi and H3 (long range) models than between any of them 
and the H2 model. This is supported by comparison of the spatial structures (snapshots) obtained at low temperatures after 
some 4 x 10^ MCS (see Figure |2| • We see that for H3 and Hi models a simple structure of the Neel type antiferromagnet is 
present, while for the H2 model each particle has vacancies as its NN and particles of the opposite charge as its NNN.. We see 
also that for short range interactions more particles are aggregated into a single cluster. Development of the spatial structure 
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FIG. 3: Time dependence of the average number of clusters in the three models. Left upper panel - low temperatures, right upper panel - 
medium ones and bottom panel - high temperatures. 
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FIG. 4: Time dependence of the average cluster size for the H3 model at low temperature. 

and its dynamics is shown in Fig|3]At low temperatures the patterns for the H3 and Hi or H2 models are different. In the H2 
model a single cluster is rather fast formed; it may break into two smaller ones, but they merge rather fast to form again one 
cluster. For the Hi model formation of a single cluster takes a longer time, but then the behavior is similar to the H2 model. 
We may therefore state that the asymptotic state for the two short-ranged models {Hi and H2) is a single cluster. Formation of 
a single cluster for the H3 model, as seen from Fig|3] is more difficult. We have run simulations in that case till 3 x 10^ MCS, 
and as shows Fig|4]even after such long time the single cluster structure is nor very stable. The reason is that with short ranged 
interactions particles which do not belong to a cluster are basically free. They move randomly and when they come close to a 
cluster boundary they quickly join the cluster and stay there. Particles in the H3 model are never really free and they experience, 
often contradictory, forces coming from many neighbors. Single particles in the H-s model are rare since they quickly form pairs 
or quadruples, which move through the system at a slower rate. 

At medium temperatures only interactions to NNN could lead to formation of a large clusters, which are nevertheless very 
unstable and quickly dissolve into smaller ones. In the case of the H3 and Hi models there is always a large number of rather 
small clusters and the respective curves lie very close one to another. Finally for high temperatures, well above Tc, as expected, 
there is no order whatsoever and no differences among the models. The dependence of the average energy on time is presented 
in Figures |5] At low temperatures we see that the energies of all models are quite different. The lowest is the energy of the 
Hi model, slightly above it lies the energy H2, while the energy in the H3 model is the highest. Most of the particles at such 
temperatures, belong to a rather well organized cluster, or clusters, with the lowest possible energy. Particles from the boundaries 
of the clusters do not have all partners, hence their contribution to the total (negative) energy is smaller. For the Hi model, in 
general, only one or two partners are missing. For the H2 model a particle interacts with more partners, hence the loss for the 
boundary particles is heavier and the total energy is raised more than for the Hi model. The same, but to a larger extent, is 
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FIG. 5: Time dependence of the average energy in the three models. Left upper panel - low temperatures, right upper panel - medium ones and 
bottom panel - high temperatures. 



true for the Hj, model. For higher (medium) temperatures all energies are shifted, the order is preserved, however the curves 
for Hi and H2 models are closer together than to the curve for the Hj, model. From Figl3]we see that at medium temperatures 
large clusters are observed only for the H2 model, hence the average energy in this case is less affected. Similarly for the Hy, 
model, where the change from the clusters size at low and medium temperatures is not dramatic. Raising the temperature has the 
largest effect on the Hi model, where the large (single) cluster dissolved at higher temperatures into many smaller ones. There 
are many more particles on the borders of the clusters and the raise in the average energy is significant. Finally, at temperatures 
well above Tc, all curves lie close together. As seen from Fig. [3] at those temperatures all larger clusters disappeared and only 
conglomerates of several particles remain. 

Mean square displacement as a function of time is shown in Fig|6]on a log-log scale for the three models and three sets of the 
reduced temperatures. It has always, at least after some initial period, a power type character: 

(r2(t))~r. (6) 

At low temperatures a = 1.0 for the Hi and H2 models, hence diffusion has normal character, while for the Hj, model we have 
a = 0.60, indicating subdiffusion. The reason for this difference lies in contradictory forces acting on separate particles in the 
H^ model, as well as in the fast formation of doublets and quadruplets. As we have shown before fl2l,T3] also in the Hi and 
H2 models subdiffusion is observed, but at temperatures lower than those reported here. For intermediate and high temperatures 
and after some IC* MCS diffusion has normal character in each case. 

From the msd data the AiThenius plots, i.e. the dependence of D* on Ji/T could be constructed using Eq.(|4ll. In Fig|2] 
on left panel we present D* as a function of the reduced temperature r. As could be expected from the previous data and its 
analysis, diffusion in the H^ model is slower, except at very high temperatures, where the particles could be treated as free. The 
dependence of D* on the temperature looks however different if we return to the unsealed temperature. Corresponding plot of 
D* versus temperature (divided by the interaction constant Ji) is shown in Fig|2]on right panel. Since the critical temperature 
for the Hj, model is much lower than the ones for the Hi and H2 models, at a given unsealed temperature, diffusion of particles 
subject to long ranged interactions is much faster 

On the last plot (Fig. (H))) we present comparison of the two-point correlation function G, measuring correlations between 
particles which are NN, and defined as 

^ N & 
i=l 3 = 1 

The same correlation function could be calculated analytically, using the Oguchi method |fl6t], where the contribution from 
interactions of two NN particles is calculated exactly, while with all others is replaced by a mean field. Derivation of the formula 
for G in the Oguchi method is described in il2tl . Here we present only the final formulae 

dlnZ sinh K 

~ dK ^ YT2^^Tx^~coshK ■ 
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FIG. 7: Arrhenius plot of the tracer diffusion coefficient on scaled temperature (left panel) and unsealed (right panel) for the three models. 



where Z is the partition function and x is 



1 - 2c 



2(1 -c) coshiC 



4c(l - c) coshi^ 



(l-2c)2 



(9) 



K — Ji/kT. c is the coverage, which here equals 0.2. The curves presented in Figl8]show similar character, except for low 
temperatures, where mean-field type theories are known to be not very reliable since they miss NN interactions which give 




FIG. 8; Comparison of the two-point correlation function obtained for the Hz and Hi models via simulations and calculated analytically 



7 



significant contribution to correlations at low temperatures. 

IV. CONCLUSIONS 

We have discussed the role of the range of interactions on static (spatial patterns, order parameter) and dynamic (msd, average 
cluster size, average energy) properties for the model introduced by us earlier llT2llT3ll . We consider here three versions of the 
model. In the first one (Hi) interactions are restricted to NN, in the second one (i?2) they are extended to NNN. In the third one 
(Hs), a particle interacts with all other particles not farther from it than half the linear size of the system, which with periodic 
boundary conditions used here means that one particle interacts practically with all other. The strength of the interactions falls 
off with the distance between the particles. We have shown that the static properties of the H3 and Hi models do not really 
differ - both exhibit the same kind of "antiferromagnetic" ordering at low temperatures. Order parameter dependence on the 
temperature scaled by the respective critical one shows the same pattern for the three models. Similarity of the two point (to 
NN) correlation function behavior indicates that the long ranged interactions are, specially at low temperatures, screened. 
Situation is however quantitatively different for dynamic features. Here taking a short (to NN or NNN) or long range range 
interactions yields different results. Organization of the system into one or two clusters takes much longer time for the H3 
model. Moreover such clusters are less stable, since there are more "unsatisfied" interactions for the border particles. Transition 
temperature from the ordered to disordered phase is much smaller for the than for the short range models, which also shows 
that the clusters formed in the H3 model are less stable. Since Tc for the H3 model is much lower, at temperatures (unsealed) 
which are still below Tc for the Hi and H2 models, particles with long range range interactions are nearly free. Hence their 
mobility is larger than those in the Hi and H2 models, still forming clusters. At low (scaled) temperatures diffusion is slower 
in the case of the H3 model, as practically on each particle act forces from many other particles and there are few, if any, free 
particles which could realize ordinary diffusive motion, characterized by simple {r^{t)) ^ t dependence. Instead, subdiffusion 
with an exponent equal 0.60 is observed. A system with long range range interactions has a tendency to organize itself, even at 
low (scaled) temperatures, in a set of smaller clusters than in the Hi or H2 models. Mobility of such groups of particles is lower 
than of separate particles encountered in the systems with short range interactions. Therefore, although finally, like in the other 
models, a single cluster appears, it takes, for the system of the size studied here, at least 10 times longer. 

Our study does not necessarily mean that when modeling real systems of charged particles diffusing on a surface, always infinite 
range interactions should be taken. First of all, static properties are reasonably well described by short range models. Moreover, 
as we have demonstrated, the differences between the three versions of our model, have generally quantitative rather than 
qualitative character. Therefore screening of the Coulomb like potential has, in the considered model, an important influence, 
reducing in some cases significantly the effective range of interactions. 
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